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The Child-Langmuir law is investigated numerically using a fully electromagnetic particle-in-cell 
code for a closed axially symmetric diode. It is found that the average current transmitted to the 
anode (Jtr) increases with the injected current (Jin) even after the formation of virtual cathode in 
both the non-relativistic and relativistic cases. The increase is found to be a power law, Jtr ~ JiV^- 
In other words, the time averaged fraction / of electrons reaching the anode varies with the input 



current as, / ~ Jj^f 
possibility of asymptotic saturation is also discussed 



Jjjj where /3 < 1. In contrast, for an infinite parallel plate diode, / 



The 



I. INTRODUCTION 

The Child-Langmuir law is the cornerstone of electron 
flow in diodes. It gives the maximum current density, 
>^Cl(1)7 that can be transported from an infinite planar 
cold cathode (velocity of emission is zero) at zero poten- 
tial to an infinite planar anode parallel to the cathode 
at a distance D and potential V. In the non-relativistic 
case, this is given by [1,2] 



^cl(I) 



V^J_ /2e 
9^ \m 



1/2 



(1) 



Here, e is the magnitude of the charge and mo the rest 
mass of an electron. 

The Child-Langmuir law is a result of the space-charge 
effect. As electrons are emitted from the cathode, they 
gain velocity from the imposed field. At the same time, 
they experience a repulsive force due to the presence 
of other electrons. However, all the injected current is 
transmitted to the anode as long as the injected current, 
Jin, is less than a particular value, called the critical cur- 
rent JcR- When Jjn exceeds Jcr, the net force is such 
that some of the electrons, instead of moving towards 
the anode, start moving backwards as if reflected from a 
"virtual cathode" . As a result, only a part of the current, 
Jtr, is transmitted. The position of the virtual cathode 
depends on Jin and on the initial velocity, vq, of the emit- 
ted electrons. In one dimension, Jcr — Jch{^) in case 
of cold emission i.e. in case vq = 0. Also, the virtual 
cathode is formed on the cathode as soon as the injected 
current exceeds Jcr and remains there even when Jin 
is increased. The transmitted current Jtr in this case 
is the Child-Langmuir current Jcl(1) for any value of 
the injected current above Jcr. The transition in the 
transported current is therefore sharp: Jtr = Jin if 

Jin < Jcr, Jtr — Jcr if >^in > ^cr- 

The flow of current below the critical value is steady 
i.e. time-independent, but is oscillatory above it. The 
value of the critical current Jcr is thus that value above 
which the system does not have a steady state. This 
criterion may be used to evaluate exact analytical ex- 
pression for Jcr in one-dimensional planar geometry in 



case the kinetic energy of the electrons is non-relativistic. 
The analytical expression for Jcr may also be derived in 
the ultra-relativistic limit. 

Since, as mentioned above, the system does not ap- 
proach any steady state if Jin > Jcr, its behaviour for 
the input currents exceeding the critical current is de- 
scribed by time-dependent equations. However, it is not 
possible to solve the time-dependent equations analyt- 
ically exactly even in one-dimensional planar geometry 
[3]. Solving time-dependent equations is avoided in the 
so called classical theory which, for Jin > Jcr, assumes 
a steady state which takes account of reflections from 
the virtual cathode. The one-dimensional planar model 
is analytically exactly solvable within the framework of 
the classical theory. In particular, it leads to the Child- 
Langmuir law. Comparison with numerical work shows 
that the predictions of the classical theory are close to 
the time- averaged values of physical quantities. 

Over the past decade, several studies have been carried 
out to extend the Child-Langmuir law to two and higher 
dimensions [4,5]. The motivation for this is twofold. In 
the first place, this helps to benchmark particle-in-cell 
(PIC) simulation codes. More importantly, the existence 
of a space charge limited saturation current implies that 
one need not bother about the material dependent cath- 
ode characteristics (the Richardson-Dushman law in case 
of thermionic emission or the Fowler-Nordheim law [7] for 
field emission) so long as the injected current is more than 
what gives rise to the limiting value of the transmitted 
current. This condition considerably simplifies the simu- 
lation of devices that use a diode. 

Most studies in two-dimensions have centred around 
cold emission from a finite emission area in an other- 
wise infinite parallel plate geometry [4-6] . Note that the 
Child-Langmuir current is defined in these studies as the 
current at which oscillatory behaviour sets in (the crit- 
ical current, Jcr ). Both numerical and analytical re- 
sults indicate that this critical current in two-dimensions 
(Jcl(2)) is higher than the Child-Langmuir current in 
one-dimension (Jci(l)). Further, there exists a relation 
between the two [5]: 
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Here W is width of the emission strip, D is the anode- 
cathode separation and ao — 0.31. The question of Hmit- 
ing current (the maximum current transported across the 
gap) has not however been investigated by these authors. 

Recent experiments using photoemission from a cath- 
ode have yielded some interesting results [8] . The experi- 
ments show that it is possible to transport a much higher 
current than that predicted by the one-dimensional law 
for a short emission pulse or limited emitter area and 
(more importantly) the current may be increased even 
after formation of the virtual cathode. This leads to 
a distinction between the critical current (the minimum 
current required for virtual cathode formation) and the 
limiting current. In the range of laser intensity stud- 
ied in these experiments, no limiting current (maximum 
current that can be transported for a given potential dif- 
ference and gap separation) was established although the 
authors suggest the existence of an asymptote. 

The question of limiting current in two and higher di- 
mensions thus assumes a great significance since its non- 
existence in realistic situations can be a major stumbling 
block in simulations. With this in mind, in this commu- 
nication we investigate numerically the Child-Langmuir 
law in an axially symmetric cavity. We demonstrate that, 
for axially symmetric closed geometries with a finite emit- 
ting area, (i) the current increases even after formation of 
the virtual cathode, and (ii) the time-averaged fraction 
of electrons reaching the anode varies as f ^ with 
/? < 1 thereby establishing that the transmitted current, 
<^TR = fJiN, increases as a power law beyond the critical 
current. 

The plan of the paper is as follows: We first review 
the one-dimensional non-relativistic case and, within the 
framework of the classical theory, show analytically that 
above the critical current, the fraction of charges reaching 
the anode varies inversely with the injected current (/ ^ 
Jj^)- A similar study for the relativistic case is presented 
next followed by our numerical results for infinite planar 
diode and axially symmetric cavity diode. Finally our 
results are summarized in the concluding section. As a 
matter of convenience, we shall use the notations J, v 
and p to refer to the magnitude of the current density, 
velocity and charge density respectively. 



shown that the potential 4>{z) at the plane at the position 
z in that case is the solution of the Poisson equation 



(4) 



in which the charge density p{z) is determined by the 
condition that the current density is independent of space 
so that p{z)v = Jin whereas the velocity ?; is determined 
by the energy conservation law fo = ^mo — e0(z) where 
£o = ^'tiqWq — e(/)(0) is the initial energy of an electron at 
the surface of the cathode. Under the condition of cold 
emission, Vq = 0, and due to the boundary condition (3) 
at z = 0, it follows that £q = 0. Hence, the law of energy 
conservation in the steady state reads 



£ = ^niov'^ — e(j){z) = 0. 
Equation (4) then assumes the form 



d^ 
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where 



C = 47rJiN(m/2e)i/2. 



(5) 



(6) 



(7) 



Equation (6) is to be solved under the boundary condi- 
tions given in (3). 

The steady state may alternatively be found by solving 
the equation of motion for electrons [9]. That equation, 
in the Llewellyn form [3,10], reads 



d3^ = -^Mt) 
dt'^ mo 



(8) 



where to is the time at which the electron is injected 
(initial time), mo the rest mass of an electron, Jt = 
dEo{t)/dt+4:TTJ is the total current and Eq is the electric 
field on the surface of the cathode. 
Eq. (8) can be integrated to get 
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II. CRITICAL AND LIMITING CURRENTS: 
ONE-DIMENSIONAL NON-RELATIVISTIC CASE 



In the steady state, Eo{t) = Eo{to) and dEo/dt = 0. 
Also, using the fact that below the critical current, no 
electron is reflected back, we have J = Jin. Thus, 



A. Determination of Critical Current 

Consider two infinite plates placed at z = and z = D 
at fixed potentials given by 



(l){z = 0) = 0, (l){z = D) = V. 



(3) 



Assume that the system evolves to a steady state. In that 

state, the current at every plane between the plates ap- 
proaches a value which is independent of time. It can be 



z{t, to) = v{t - to) - - — Eoito){t - tof 
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and 



v{t, to) =vo- —Eo{to){t - to) + l^±^(t - to)^. 
mo 2mo 



(11) 
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For the case of cold emission {vq = 0), (5) holds giving 
V = (2e(/)/mo)^/^ so that the final velocity at the anode is 
Vd = (2ey/mo)^/^, where (/)(D) = V. Let T denote the 
time that an electron takes to transit the cathode-anode 
distance {D). Then, with z = D, v = vd, t — to = T, 
Eqs. (10) and (11) assume the form 
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These equations determine the two unknowns, Eo and 
T. On eliminating Eq, the equation for T in terms of the 
scaled transit time T is 



g 

T — aT + a 







(14) 



where a = 27V^/^/4P, T = T/Tq, To = D/{2mo/eVy/^ 
and /3 = 97rJiN£'^(mo/2e)i/2. Note that To is the transit 
time in the absence of space charge. The condition that 
all the roots of (14) be real may be shown to be a > 27/4. 
For a < 27/4, two of the roots are complex while one 
root is real. The real root is however negative and hence 
inadmissible as a transit time. Thus a — 27/4 marks the 
critical current. On substituting the values of a and /?, 
the expression for the critical current is found to be given 

by 
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Also, note that at Jin = JcRi T = 3/2 and hence Eq = 
0. Thus, at the critical current, the electric field at the 
cathode vanishes. 



where (f)= {^ + £o/e-)/V, z = z/D, and 
97rJiN£'^ 



a = 
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Both these equations are of the form d?(l)/dz'^ = A/cp^^^ 
and can be cast in the form 



dz \ dz / dz 



The solution of this equation reads 

2 ^ 
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(20) 



Since the electrons come to a stop at z = Zm, the law of 
conservation of energy {£o = mov"^ /2 — e<f)) implies that 
(t>{zm) = 0. Also, the field on the electrons at ^; = is 
assumed to be zero. Hence 



d^ 

dz 



= 0. 



(21) 



On using the abovementioned values of 4> and d4>/dz at 



Zm, it follows from Eq. (20) that B = 0. Thus 
d^ 



dz 



(22) 



where minus sign holds in the region < z < z^ whereas 
the plus sign is for the region z^ < z < 1- That is 
because (f>{z) decreases from ^ = to ^ = 2m so that the 
electrons arc decelerated and increases beyond Zm so that 
a fraction of electrons is transmitted after they come to 
rest at z = Zm- The solution of (22) yields 



B. Determination of Limiting Current 
Theory 



Classical 



The above analysis holds good for Jin < JcR where 
the injected and transmitted currents are equal. Beyond 
the critical current, electrons are reflected and the anal- 
ysis breaks down. In order to find out what happens 
beyond the critical current, it is useful to consider yet 
another steady state model which assumes that a frac- 
tion / of the injected electrons is transmitted while the 
fraction 1 — / is reflected at some point z = Zm so that 
the transmitted current is Jtr = / Jin [2,11]. Though we 
have assumed so far that the total energy £o of an elec- 
tron entering the diode is zero, it is instructive to carry 
the classical theory for injected currents greater than the 
critical currents for the case of non-zero £o- To that end, 
on invoking the discussion following Eq.(4), the Poisson 
equation in the two regions, < z < Zm and Zm ^ z < D , 
then assumes the form 



d^ 

dz2 

d2^ 

di2 



4a(2 - /) 

9v^ 
4a/ 
9^^ 
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^^/' = lA{z-Zmr. 
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On using the appropriate values of A in the two regions, 
the solution can be expressed as 



P/' = ai2 - f){z - z^Y' 

p/"^ = af{z-Zmf, Zr, 



< Z < Zm, 
<Z<1. 



(24) 



Now, on applying the boundary conditions (3), the equa- 
tions in (24) yield 

3/2 
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-) 

y + gp/e ^'/' 

V 



, = a(2 - f)zl, 



= af{l - Zmf. 



These equations show that if, £"0 = then 

•Zm = and af = 1. 



(25) 



(26) 



The position of the virtual cathode is thus on the cath- 
ode when the initial energy of an electron is zero but is 

away from it otherwise. On substituting for a from (18) 
in (26), the transmitted current is found to be given by 
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Note that this value of the transmitted current is inde- 
pendent of the injected current and that this expression 
holds as long as the injected current exceeds the critical 
current. Eq. (27) is thus the limiting current that can 
flow through the diode. Note also that the limiting cur- 
rent (27) is identical to the critical current Jcr given in 
(15). 

It should be emphasized that the steady state model 
for Jin > Jcr is phcnomcnological. As is revealed by 
the numerical solution of the time-dependent equations 
for injected currents above Jcr, the behaviour of the 
system is oscillatory i.e. it does not approach a steady 
state. However, the steady state values derived above are 
found to be in close agreement with the time-averaged 
values of the numerical results. 



III. CRITICAL AND LIMITING CURRENTS: 
ONE-DIMENSIONAL RELATIVISTIC CASE 

Finding an explicit expression for the critical current 
in case of relativistic electronic energies is generally a 
formidable task. However, an analytical expression for 
the critical current may be derived in the ultrarclativistic 
limit in case the initial velocity of the electrons is zero 
if, following the results of the non-relativistic treatment, 
we identify the critical current for zero initial velocity 
as the one for which the electric field at the cathode is 
zero. We will see that, as in the non-relativistic limit, 
this current is the same as the limiting current predicted 
by the classical theory for relativistic energies. 



A. Critical Current 

The expression for the critical current in the relativis- 
tic case may be derived more conveniently by starting 
from the Poisson equation (4). The charge density p in 
(4) is determined by pv = Jin and v by the relativistic 
energy conservation law 



moc^i'jiz) - 1) - e(p{z) = 0, 



(28) 



where ^{z) = (1 - v^{z)/c'^) ^. Eq. (4) may then be 
rewritten as 



dz^ 



(29) 



where K = A-k Jm{e / niQC?) . The boundary conditions 
(3) in this case assume the form 



^[z = 0) = 1, ^{z = D) = l + 



eV 
moc^ 



ID. (30) 



Following the non-relativistic case, we assume that 

when the initial velocity of the electrons is zero, the criti- 
cal value of transmitted current is attained when the field 



at cathode vanishes. Hence, assuming d(j}{z)/dz = at 
z = and invoking also the boundary condition (30) at 
z = Q, the integration of (29) gives 



dz 



= 2K{Y - 1) 



2 _ 1^/2 



(31) 



On integrating this equation with the boundary condi- 
tion (30) at z = it may be shown that the current 
density (which, as explained above, is to be identified as 
the critical current density) is given by 



Jcr 



87762;' 



where 



I{l) = j\-"\l-l-')-'/'d^. 



(32) 



(33) 



Eq. (32) determines 7 implicitly as a function of z and Jqr 
in case the field at the cathode is zero. However, Jcr is 
still an unknown. It may be determined by applying the 
yet unused boundary condition B.t z = D leading to the 
relation 
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where 7d is given by (30) in terms of the applied voltage 
V . Hence, (34) determines the critical current in terms 
of known quantities. 

The function 1(7) in (33) may be evaluated as follows: 
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where F{x) = F(-|, — ^, |; x) stands for the Hypcrgco- 
metric function. Now, let 7 = 7/5 and let the voltage 
be such that eV » rrioc^ so that » 1- Hence, 



F(a,6,c;l) = 



r(c- 



6)r(c) 



r(c-a)r(c-6)' 



(36) 



and r(3/4) = 1.2254167024 [12], the expression for the 
critical current in the ultrarelativistic limit turns out to 

be given by 



-^CR 



2neD^ 
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0.8471 



(37) 



This is in agreement with the result of Jory and Triv- 
elpiece [13] derived under identical boundary conditions 
though these authors do not identify (37) as the critical 
current. 
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B. Limiting Current 



IV. NUMERICAL RESULTS 



As in the non-relativistic case, we assume that when 
Jin > JcR, the steady state is characterized by a frac- 
tion / of particles transmitted beyond the position z„i 
at which the velocity becomes zero so that the current 
in the region between the cathode and the position Zm 
consists of two parts, the injected current density Jin 
moving away from the cathode to the virtual cathode at 
Zm and the reflected part (1 — /) Jin moving to the cath- 
ode from the virtual cathode. The Poisson equation in 
the two regions can then be written as 



(2 - /) 
fK-t 



d^ 

dz2 

d^ _ 

dz^ ~ .J^f^ 



^< Z < Zn 



< z< D, 



(38) 
(39) 



where we have assumed that the energy of the injected 
electrons is zero. Energy conservation thus leads to the 

condition (7— l)moc^ — etp — O. The two equations above 
are of the same form, whose solution is 



dz 



c 



(40) 



where K = {2 - f)K for < z < Zm and K = fK for 
Zm l£ z < D. Note that, at z = Zm, 7=1, d7(z)/dz = 0. 
Thus C = 0. Keeping in mind that the potential de- 
creases in (0, Zm) and increases in {zm, D), the solutions 
in the two regions read: 



d7 



1 (7^-1)1/* 



{2fKf'\z - Zm), Z > Zm. (41) 



Like in the non-rclativistic case discussed in the last sec- 
tion, the analysis may be carried for non-zero values of 
initial energy Sq. Restricting, however, to the case of 
£0 = with zero initial velocity (so that j{z = 0) = 1) 
and zero potential at the cathode, the boundary condi- 
tions give Zm = and 



in 



nioc 
Sz^ire 



I\l), 



(42) 



where 1(7) is as in Eq. (33). Eq.(42) with / Jin JcR 
is the same as (32). On applying the boundary condi- 
tion 7 = 7£) at 2: = £> we get the current at the anode. 
Hence, the transmitted current /Jin is the same irre- 
spective of the value of the injected current Jin as long 
as Jin is above the critical value at which virtual cathode 
is formed. The critical value of the current is, therefore, 
the same as the limiting current. Note from (42) that the 
fraction of electrons transmitted varies as / ~ Jin^- 



The 1-D electrostatic analysis reviewed above indi- 
cates that, for cold emission, the critical and limiting 
currents are identical in both non-relativistic and rela- 
tivistic cases. The chief motivation of this study is to 
investigate whether this holds in higher dimension. To 
this end, we shall use the fully electromagnetic particle- 
in-cell code, SPIFFE [14], and limit ourselves to axially 
symmetric diodes. For completeness and comparison, we 
shall first present our numerical results for 1-dimension. 



A. One-Dimension: Sheet Model 

In practical terms, the 1-dimensional analysis applies 
when the cathode and anode are parallel plates of di- 
mensions much larger than the separation between them 
with uniform emission from the surface of the cathode. 
A description, convenient for numerical calculations, is 
to discretize continuous electronic fluid as sheets of uni- 
form charge density parallel to the surface of the two 
electrodes [3]. The sheets move in the direction perpen- 
dicular to their surface. The position Zk of the k*^ sheet 
at time t is governed by the equation [16] 



dt 



2 Qfrel 



M 



(43) 



This equation takes into account the electrostatic re- 
pulsion between the sheets. Here Zk = Zk/D is the 
scaled position of the k^^ sheet, (3k = vt/c, aroi = 
47r Jin J'^(9e/2moc'^)(c/tio)"^, vq is the initial velocity, N 
is the number of sheets launched per unit transit time in 
the absence of space charge and M is the total number 
of sheets present in the diode at any instant of time. 

For the numerical calculations presented here, the in- 
tegration time step At is D/(200vo), D = 0.8 cm while 
N = 2000 for non-relativistic calculations {V = 250 kV) 
and N = 10000 for rclativistic calculations {V = 2 MV). 
These parameters were chosen to satisfy convergence re- 
quirements. 

As stated before, the system approaches a steady state 
if the injected current is below its critical value whereas 
the solution is oscillatory for currents above that value. 
Since the classical theory is based on the assumption of a 
steady state even above the critical current, we compare 
it with the time-averaged behaviour of the physical quan- 
tities if the injected current is above the critical current 
[17]. 

Fig. 1 is a plot of the transmitted current averaged over 
10ns as a function of the injected current for V = 250 
kV. It is clear that Jtr attains saturation at Jin = 
4.3 MA/m^. In order to quantify the saturation, we have 
studied the (time averaged) fraction / of electrons reach- 
ing the anode as a function of the injected current. Ac- 
cording to the analysis presented in section II, / ^ J] 



Fig. 2 confirms that / ~ with /? = 0.999 



IN 



5 



4.5 r 



4 - 1 
3.5 - 1 



Jtr 3.0 



2.5 
2.0 
1.5 



40 60 
Jin 



(2 MV) leads us to an identical conclusion - the critical 
and limiting currents arc indeed identical. Fig. 3 provides 
confirmation of this. Note that the ultra-relativistic anal- 
ysis, presented in section III predicts Jqr — 80 MA/m^ 
while the observed value is 78.3 MA/m^. The saturation 
has also been studied using the averaged fraction of trans- 
mitted electrons as a function of the injected current. It 
is found that / ~ Jf^-^^^. 

Thus, in the 1-dimensional case, the critical and lim- 
iting currents are indeed identical and, above the critical 
current, the fraction of electrons transmitted varies as 



FIG. 1. Time averaged transmitted current density plotted 
as a function of the injected current density for the 1-dimcn- 
sional case. Both are measured in MA/m^. Note the sudden 
transition at 4.3 where saturation sets in. The theoretical 
value using Eq. (1) is 4.5. 
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FIG. 2. A log-log plot of the average fraction of transmit- 
ted electrons versus the injected current for the 1-dimensional 
case. The linear fit confirms that / ^ with /3 = 0.999. 
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FIG. 3. Time averaged transmitted current density plot- 
ted as a function of the injected current density for the 1-di- 
mensional relativistic case. Both are measured in MA/m^. 
The critical current density is found to be 78.3 MA/m^ after 
which saturation sets in. The theoretical value using (37) is 
80 MA/m^. 

A similar study of the 1-dimensional relativistic case 



B. Axially Symmetric Diode: PIC 

For the axially symmetric diode, the numerical calcu- 
lations were performed using the fully electromagnetic 
particle-in-cell code SPIFFE [14] . The basic algorithm is 
described below [15]. 

1. Create the computational mesh and assign metal or 
vacuum points according to the geometry specified 
in the input file. 

2. Evaluate imposed fields due to potentials applied 
on metal surfaces 

3. Inject macroparticles with specified charge and ve- 
locity using the "over-injection" method. 

4. Distribute the charge and velocity to the mesh 
points using a standard weighting scheme. 

5. Solve Maxwell's curl equations with the speci- 
fied boimdary conditions (Dirichlet/Neumann) us- 
ing the Finite Difference Time Domain (FDTD) 
method. 

6. Find forces at the position of the particles by inter- 
polating the field values at the adjacent grid points. 

7. Using these forces, find new position and momen- 
tmn of the particles. 

8. Remove particles that reach the end of the simula- 
tion region or hit a metal surface. 

9. Repeat 1 to 8. 

After a specified number of time steps, it checks if 
Poisson's equation is satisfied within a specified "error 
charge". If not, it corrects the electric field after solving 
Poisson's equation using the error charge. 

The geometry of the diode consists of a hollow cylinder 
of radius 12.5 cm with the cathode and anode plates of 
radius 12.5 cm placed on either end with a separation 
D = 0.8 cm between them. The cylinder and anode plate 
are connected and maintained at a potential, V. The 
cathode plate (on the left in fig. 4) is grounded and the 
emitting area is restricted to a radius R=3.5 cm. In each 
case, the integration time step was chosen to be 0.0002 
ns. Convergence was checked against the mesh spacings 
in r and z as well as the charge per macro-particle. 
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FIG. 4. The geometry of the diode and a typical plot of the 
particles for Jin > Jcr- Here Jin = 6MA/m^ and V = 250 
kV. 
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FIG. 7. The electric field after 2 ns. Note that the field at 
the cathode is nearly zero. 
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FIG. 5. The phase space plot corresponding to fig. 4. Note 
that pz is measured in units of moc. In the absence of space 
charge, pz = 1-1 at the anode for V = 250 kV. 
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FIG. 8. The time averaged current density close to the an- 
ode plotted against the radius. Note that the emitter radius 
was 0.035 m. 
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FIG. 6. The imposed electric field, Ez,imp due to the ap- 
plied potential 



The simulation geometry along with a typical plot of 
the particles in configuration space for Jin > Jcr is 
shown in fig. 4. The phase space plot corresponding to 
fig. 4 is shown in fig. 5. 
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FIG. 9. The locally time averaged current density as a func- 
tion of time. Note that it is steady after the initial transients. 

Figs. 6 and 7 are plots of the electric field at the 
beginning of the simulation and after 2 ns when the tran- 
sients have settled. Here Jin = 6 MA/m^ and r is fixed 
at 0.01 m. Note that in fig. 6, the field is due to the 
applied potential alone while in fig. 7 the field is a super- 
position of the applied field and the electromagnetic field 
generated by the charges. E^^ changes little after acquir- 
ing the form shown in fig. 7 characterized by its nearly 



7 



vanishing value at the cathode. 

The radial profile of the current density for Jin = 
6 MA/m^ is shown in fig. 8. The plotted value of the 
z component of the transmitted current density, J^, is 
averaged over 10 ns and the section is taken close to the 
anode. The variation of the radially averaged current 
density with time is shown in fig. 9. In order to reduce 
fluctuations, a local time averaging has also been per- 
formed. 
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current (Fig. 11). The power law behaviour is evident 
with / ~ -^in''^^- Thus the transmitted current increases 
beyond the critical current as Jtr = / Jin ^ -^in*^*- 



Jtr 40 




20 40 60 80 100 120 140 
Jin 



FIG. 12. Time averaged transmitted current density plot- 
ted against the injected current density for F = 2 MV. Both 
are measured in MA/m^. 



FIG. 10. Time averaged transmitted current density plot- 
ted against the injected current density for the axially sym- 
metric diode with V = 250 kV. Both are measured in MA/m^. 
Unlike the 1-dimensional case, no saturation can be observed 
for Jin > Jcr- 
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FIG. 11. A log-log plot of the average fraction of trans- 
mitted electrons versus the injected current for the axially 
symmetric diode with V = 250 kV. The linear fit for Jin > 5 
confirms that / ~ Jjn'' with (3 = 0.92. 

We now present the central result of this communica- 
tion. For the non-relativistic studies, the potential con- 
sidered was 250 kV for which Jcl(1) = 4.57 MA/m^ 
and Jcl(2) = 1.0354 Jcl(1) = 4.73 MA/m^. The in- 
jected current was varied in the range < Jin < 100. 
While criticality sets in at around Jcl(1), there appears 
to be no saturation in this range (sec Fig. 10). While 
Jtr — Jin for Jjn < Jcr-, the transmitted current in- 
creases as a power law beyond the critical current. This 
is evident from a plot of the average fraction of elec- 
trons reaching the anode as a function of the injected 
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FIG. 13. A log-log plot of the average fraction of transmit- 
ted electrons versus the injected current for V=2 MV. The fit 
is made for Jin > 50 and the value of /J is 0.873. 

A similar study for a relativistic case (2 MV) leads 
us to an identical conclusion - the transmitted current 
increases beyond the critical current as evident from 
Fig. 12. The power law behaviour can be seen from a 
plot of ln(/) vs In(JiN) (Fig. 13) and for this particular 
case, /3 = 0.873. Thus, Jtr ~ JiV^^- Note that in the 
relativistic case, there is a strong pinching effect due to 
the self magnetic field. The observed value of the critical 
current, Jcr, is 47 MA/m^ which is much less than the 
1-dimensional prediction (80 MA/m^). 

These results do not preclude the existence of asymp- 
totic saturation. Rather, a fit of Jtr using the form 
•^SAT — ^/-^iN provides a marginally better fit than 
the power law. However, a unique set of values for 
{JsAT, A,"/) is hard to determine with the existing data 
as several sets exist with nearly the same error (sum of 
squares). Nevertheless, it is clear from these studies that 
Jtr does not saturate for finite values of injected currents 
beyond the critical value. 
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V. CONCLUSIONS 

The existence of a limiting current in an axially sym- 
metric cavity diode has been investigated. It is found 
that Jtr does not saturate with an increase in Jin be- 
yond the critical value characterised by the onset of re- 
flection. This is unlike the case of one-dimensional diode 
which exhibits saturation of Jtr as soon as Jin exceeds 
JcR. Our main result is the power law behaviour, es- 
tablished through numerical computations between the 
transmitted current Jtr and finite values of the in- 
jected current Jin in an axially symmetric diode oper- 
ating above the critical current Jcr. The exponent de- 
pends on the voltage. We have also explored the pos- 
sibility of asymptotic saturation and found that there 
exist sets of parameters (Jsat,^, 7) for which the form 
Jtr — JsAT"^/ Jin provides a marginally better fit than 
the power law Jtr ^ Jj'^. 

The findings here thus have an important bearing on 
numerical simulations which generally assume saturation 
of transmitted current even in two-dimensional diodes at 
a finite value of the injected currenbt. Analytical and nu- 
merical studies of a one-dimensional diode have also been 
presented for the sake of comparison and completeness. 
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